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ABSTRACT 

We use high-resolution time-dependent numerical simulations of accretion discs around 
white dwarfs to study the structure and properties of transition fronts in the context 
of the thermal-viscous disc instability model. The thermal structure of cooling and 
heating fronts is dominated by radiative cooling and viscous heating, respectively, 
except in a very narrow precursor region in heating fronts where advection and radial 
transport of energy dominate. Cooling fronts are much broader than heating fronts, 
but the widths of both types of fronts scale with the local vertical scale height of the 
disc. We confirm that during a fair fraction of the propagation time of a cooling front, 
the structure of the inner disc is close to self-similar. The speed of heating fronts 
is ~ a few km s _1 , while the speed of cooling fronts is ~ a fraction of a km s _1 . 
We show that direct measurements of the speed of transition fronts probably cannot 
discriminate between various prescriptions proposed for the viscosity parameter a. 
A natural prediction of the disc instability model is that fronts decelerate as they 
propagate in the disc, independent of the prescription for a. Observation of this effect 
would confirm that dwarf nova outbursts are driven by the thermal-viscous instability. 
Most of our results also apply to low mass X-ray binaries in which the accreting object 
is a neutron star or a black hole. 
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1 INTRODUCTION 

It is widely accepted that the disk instability model (DIM) 
provides the correct description of dwarf nova outbursts 
(e.g Bath & Pringle 1982, Smak 1984, Lin, Papaloizou & 
Faulkner 1985, Cannizzo 1993a) and probably of soft X-ray 
transient events (van Paradijs & Verbunt 1984, Mineshige 
& Wheeler 1989, Cannizzo 1998a). In this model, the vari- 
ability of a steadily fed thin accretion disc is due to the 
occurrence of a thermal-viscous instability in the disc. For 
a certain range of mass transfer rates from the secondary, 
Mt, the effective temperature of a disc annulus lies within 
the unstable range 5000 — 8000 K, which corresponds to hy- 
drogen recombination inside the disc. The unstable annulus 
then experiences a limit cycle during which matter is pro- 
cessed at rates larger or smaller than the mean rate Mt (see 
Cannizzo 1993b or Osaki 1996 for reviews). 

Transition fronts are essential to the DIM because they 
are the link between the local instability of an annulus and 
the global evolution of the disc. The time required by a 



heating front to cross the disc is related (depending on wave- 
length) to the rise time of an outburst. Similarly, the time 
required for a cooling front to cross the disc defines the decay 
time of an outburst in the DIM (if the entire disc is brought 
to the hot state, the decay time is increased by a "viscous 
plateau" phase, see Cannizzo 1993a). Observations of the 
rise, decay and recurrence times therefore allow a direct test 
of the predictions of the DIM. 

Several attempts to determine the structure and prop- 
erties of transition fronts exist in the literature. Analytical 
studies lead only to qualitative results or to complex results 
that require numerical simulations for calibration (Meyer 
1984, Lin et al. 1985, Vishniac & Wheeler 1996, Vishniac 
1997). On the other hand, numerical studies have difficulties 
in correctly handling transition fronts because the fronts are 
very narrow and their propagation involves time and length 
scales which vary by several orders of magnitude in a single 
outburst. As discussed by Hameury et al. (1998, hereafter 
HMDLH), heating fronts are unresolved in most numerical 
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simulat ons while cooling fronts are often barely resolved, 



in particular close to the disc inner edge (see e.g. Cannizzo, 
Chen & Livio 1995). The poor resolution of transition fronts 
not only affects localized regions of the disc but it has also 
consequences for the global evolution of the disc (Cannizzo 
1993a, HMDLH). 

HMDLH have recently developed an implicit time- 
dependent numerical code which solves the disc equations 
on an adaptive grid. The code allows the resolution of very 
narrow structures in the disc at a relatively low computa- 
tional cost. Here we use this code to study the detailed struc- 
ture and the propagation velocities of heating and cooling 
fronts. We show that the width of any transition front is 
proportional to the disc pressure scale he ight i f, and not to 
(HR) 1 / 2 as proposed by Cannizzo et al. (1995). We confirm 
the existence of a self-similar regime of the hot inner disc 
during the propagation of a cooling front, and we discuss 
the possible reflections of cooling (resp. heating) fronts into 
heating (resp. cooling) fronts. 



2 DISC INSTABILITY MODEL AND 
TRANSITION FRONTS 

2.1 Disc equations 

The time evolution of an accretion disc is driven by thermal- 
viscous processes. For a Keplerian disc the viscous equation, 
describing ma ss and angular momentum con servation, can 
be written as ( Lightman, 1974[ |Pringle, 1981 ): 



9E 
dt 



3__9_ 

RdR 



R 



1/2 



OR 



(vXR 1/2 ) 



(1) 



where v is the kinematic viscosity, E is the surface density of 
the disc, and R is the distance from the central object. The 
viscosity is parameterized according to v = ac s H, where c a 
is the sound speed in the disc midplane, H is the vertical 
scale height of the disc and the parameter a (< 1) describes 
our ignorance of the viscous processes. Shakura and Sunyaev 
(1973) first proposed this prescription with a constant, but 
it was later shown by Smak (1984) that a must be larger for 
the hot ionized disc than for the cold neutral disc in order 
to obtain large amplitude outbursts in the DIM. 

In the following, we use two parameterizations often 
adopted in the literature. First, we adopt an a no t — Qcoid 
prescription and an arbitrary functional dependence of a 
with the disc central temperature T c which is similar to that 
used in HMDLH, i.e. 



log (a) = log(« co i d ) + [log(a ho t) - log(acoid)] 



1 + 



2.5 x 10 4 K 



(2) 



where Qhot and a co id are the two constant values of a on the 
hot and the cold stable branches of the thermal equilibrium 
curves (or S-curves), respectively. 

The second parameterization we consider is 



a = a {H/Rf 



(3) 



where ao is a constant (ao = 50 here). Such a depen- 
dence has been strongly advocated by several authors to 
explai n the observed exponential decays of soft X-ray tran- 
sients (Cannizzo, Chen fc Livio, 1995; Vishniac fc Wheeler 



199f). This interpretation may face some difficulties, how- 
ever (§4.6). Note that the viscosity parameter a tx R 3 ^ in 
this case (e.g. Mineshige & Wheeler 1989). 

The importance of allowing the disc outer edge, Rout, 
to vary with time in order to reproduce correct outburst cy- 
cles was emphasized by HMDLH and before by Smak (1984) 
and Ichikawa and Osaki (1992). Transition fronts are, how- 
ever, localized structures and should not be influenced by 
the physical conditions at the disc outer edge. This is con- 
firmed by a comparison between models differing only in 
whether or not -Rout varies with time, which shows that the 
structure and the properties of the fronts are identical in the 
two classes of models (except when a front reaches -Rout). In 
the following, we assume that R ou t is fixed. 

The thermal equation (describing energy conservation; 
Faulkner, Lin & Papaloizou 1983, HMDLH) is given by: 

dTo _ 2(Q+~Q- + J) _ gtT c 1 d(RV T ) 

dt 



CpT, 



HC P R OR 



where Q + and Q~ are the vertically-averaged rates of vis- 
cous heating and radiative cooling, respectively, 5R is the 
perfect gas constant, Cp is the specific heat at constant 
pressure of the gas, fi its mean molecular weight and V T 
its radial velocity (positive outward). The second term on 
the r.h.s. of Eq. (Q) represents the contribution from pres- 
sure work and will be referred to as Q p dv The third term 
(Qadv) represents the contribution from energy advection. 
The term J (or equivalently Qj) accounts for the radial 
transport of energy (viscous or radiative). The (viscous) flux 
carried by turbulent eddies with a characteristic velocity v e 
and size l c is obtained dimensionally in the framework of the 
a-prescription: 



-F c = Cp'Evc^^-lc 
OR 



2" CpE dR> 



and J is then given by: 
1 d 



J 



RdR 



(RFC 



(•») 



(6) 



A similar e xpression is obta ined for the radiative transport 
of energy (Lin et al., 1985). Other prescriptions were pro- 
posed for J but they give results which are similar, in terms 
of outburst cycles, to those obtained with Eqs. (g) and Q 
(HMDLH; see also §6). 

Simulations show that most of the disc remains close 
to thermal equilibrium during its time evolution. In that 
case, the thermal equation basically reduces to Q + — Q~ . 
Transition fronts correspond, however, to regions far from 
thermal equilibrium in which the other terms in Eq. (^) are 
no longer negligible. 

2.2 Numerical models 

We use the numerical code developed by HMDLH with 
800 grid points. This corresponds to about 100 points in 
the transition fronts, and is more than sufficient to avoid 
resolution-limited results. 

All the models predict strictly periodic sequences of out- 
bursts with (eventually) various amplitudes (see e.g. Fig. 8 
of HMDLH). The sequences differ from one model to an- 
other (in number and amplitude of the outbursts). In the 
following, we refer to such a sequence as the outburst cycle 
of a model. 
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The disc outer radius is fixed to R ou t = 4 x 10 cm. 
The disc inner radius, Ri n = 5 x 10 s cm, is equal to the white 
dwarf radius of mass Mi = 1.2 M@. A comparison of models 
with various mass transfer rates (Mt) shows that, although 
the outburst cycles differ from one model to another, the 
intrinsic structure and properties of transition fronts do not 
depend on M T . A value M T = 10" 9 M Q yr" 1 = 6.66 x 10 16 
g s _1 is used in all models presented here. 

We have considered four models to investigate the in- 
fluence of the viscosity prescription on the structure and the 
properties of transition fronts. In three of them, we use an 
«hot — «coid prescription: a hot =0.1 and a co id = 0.02 (model 
hO.l.cO.02), a hot = 0.2 and q coM = 0.02 (model h0.2.c0.02) 
and finally ahot = 0.1 and a co id = 0.01 (model hO.l.cO.01). 
The fourth model uses a = 50(H/R) 15 (model 50hrl.5). 

Global stability considerations (e.g. Cannizzo 1993b) 
and numerical simulations show that cooling fronts invari- 
ably appear in the outer disc. On the contrary, outbursts 
may be triggered anywhere in the disc, and both inside-out 
and outside-in heating fronts are observed in simulations 
(Smak 1984, HMDLH). All outbursts are of the inside-out 
type in the models presented here, except the largest ampli- 
tude outbursts of model hO.l.cO.01 and all the outburst of 
model 50hrl.5 which are of the outside-in type. 

2.3 Transition fronts 

We define the location of a front, -Rfront, as the radius at 
which the central temperature T c = 1.8 x 10 4 K in the disc. 
This value of T c is a signature of the presence of a transition 
front because it lies in the range of central temperatures for 
which an annulus is thermally and viscously unstable. The 
speed of a front, Vf ron t, is obtained by numerical differenti- 
ation Of -Rfront. 

In the following, the structure of transition fronts is 
shown each time for a specific value of Rf r0 nt- Simulations 
show that the structures of transition fronts are qualitatively 
independent of their location in the disc. (In addition, cool- 
ing fronts evolve in a nearly self-similar way.) The general 
properties of the fronts do not depend on the viscosity pre- 
scription either. 

From now on, the various terms Q/Q + in Eq. (^) used 
to describe the thermal structure of a transition front are 
defined as positive if they correspond to a local heating and 
negative if they correspond to a local cooling in the disc (i.e. 
contrary to the positive convention for Q~ in Eq. (Q)). 
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Figure 1. Structure of an inside-out heating front located at 
Rfront ~ 3 X 10 9 cm (model hO. I.c0.02). The structure is domi- 
nated by a strong outflow of gas and by viscous heating, except 
in a narrow precursor region where radial diffusion and advection 
of energy dominate. 

inside-out heating front is, by far, the viscous heating term 
Q + . Figure [l] shows, however, the existence of a substruc- 
ture in the front: a precursor region (narrower than the front 
itself) in which the energetics of the gas is dominated by 
advection and radial transport of energy. The contribution 
from pressure work is small but not negligible in the precur- 
sor region. The thermal structure of the front can be under- 
stood as follows. The hot gas which penetrates in colder re- 
gions of the disc carries a substantial amount of heat (Q a dv). 
The gas is slowed down where the transport of angular mo- 
mentum is reduced and it contracts (Q p dv)- The precursor 
region is the region of strongest gradient of T c (hence the 
large value of Qj). 

Note that the region where heating (Q + ) dominates 
over cooling is broader than the spike in E and the pre- 
cursor region corresponds to the rising part of the spike of 
E. 



3 HEATING FRONTS 

3.1 Structure of inside-out heating fronts 

Figure |l| shows the structure of an inside-out heating front 
located at R lro nt ~ 3 x 10 9 cm in model h0.1.c0.02. It is 
characterized by a strong outflow of gas which is coincident 
with a sharp spike in the E profile and a steep gradient of 
T c (or equivalently viscosity). The outflow and the forma- 
tion of the spike in E were explained by Lin et al. (1985); 
they appear because the gas which carries excess angular 
momentum from the inner disc reaches colder regions where 
transport of angular momentum is reduced. 

The dominant term in the thermal equation inside an 



3.2 Structure of outside-in heating fronts 

Figure |^ shows the structure of an outside-in heating front 
located at Rtmnt ~ 1-5 x 10 9 cm in model 50hrl.5. The 
structure is expected to be different from that of an inside- 
out heating front since the matter crossing the transition 
region no longer carries excess angular momentum from the 
inner disc. Instead, this matter now corresponds to the bulk 
of accretion in the disc, while exce ss angular mom entum is 
being freely carried away outward (Lin et al., 1985). 



The accretion of hot gas is revealed by the strong in- 
flow and the spike of E where the gas enters colder regions. 
Note the very steep gradient of T c at the interface between 
the hot and cold regions of the disc and the presence of a 
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Figure 2. Structure of an outside-in heating front located at 
Vfront ~ 1-53 X 10 9 cm (model 50hrl.5). The structure is domi- 
nated by a strong inflow of gas and by viscous heating, except in a 
very narrow precursor region where radial diffusion and advection 
of heat dominate. 

substructure in the profiles of E, V r and M (corresponding 
to the precursor region). 

The thermal structure of an outside-in heating front is 
similar to that of an inside-out heating front, despite their 
opposite direction of propagation. Viscous heating (Q + ) is 
dominant, except in the precursor region where advection 
and radial diffusion of energy dominate the energetics of the 
gas. Once again, the pressure work term (Q p dv) is small but 
not negligible in the precursor region. The precursor region 
is, however, much narrower (as compared to the total width 
of the front) and its structure more complex than in an 
inside-out heating front. 

3.3 Heating front velocity 

Figure ^ shows the speed of heating fronts for the four mod- 
els with various a-prescriptions. Each dot represents the 
speed of the front, Vf ront , at a specific location Rf Iont in 
the disc. Successive outbursts are triggered at nearly the 
same "ignition" radius during an outburst cycle in all the 
models, except in model hO.l.cO.01 for which the outburst 
cycle is made of one outside-in (not shown here for clarity) 
and several inside-out outbursts. The simulations show that 
the front velocities strongly depend on the surface density 
profile at the onset of the thermal- viscous instability. Since 
mass piles up and diffuses from the outer regions of the disc 
during an outburst cycle, the front velocities vary from one 
outburst to another during the cycle, as seen from Fig. |j| 
For comparison, the sound speed is ~ 15 km s _1 for the 
characteristic temperature T c — 18, 000 K in the transition 
front. 

The velocity profiles of outside-in heating fronts (model 



Figure 3. The four panels show the speed of successive heat- 
ing fronts in our four models with various a-prescriptions (dots). 
Outside-in (model BOhrl.S) and inside-out (three other models) 
heating fronts slow down during their propagation in the disc. 
The speed of heating fronts significantly depends on the value of 

«hot- 

50hrl.5) are very different from those of inside-out heating 
fronts (three other models). (This is not just due to the 
different a-prescriptions.) A common characteristic of the 
two types of fronts is, however, a deceleration as the front 
propagates in the disc (a characteristic of cooling fronts as 
well; see §4.2). 

It is not possible to compare quantitatively the speed of 
inside-out and outside-in heating fronts from Fig. ^ because 
a different a-prescription is used in the models where the two 
types of fronts appear. Fortunately, we obtain both outside- 
in and inside-out outbursts in our model hO.l.cO.01. The 
maximum speeds of the two types of fronts in this model 
are similar (~ 3 km s" 1 ) but outside-in heating fronts cross 
the disc faster than inside-out heating fronts because their 
speed is high over a larger range of radii. This is in agreement 
with the fact that outside-in outbursts are known to produce 
more asymmetric lightcurves (i.e. shorter rise times) than 
inside-out outbursts (e.g. Smak 1984, HMDLH) . 

The speed of heating fronts slightly increases when a co id 
decreases. Increasing the value of a^ ot results in a substan- 
tial, approximately linear, increase of Vfront which is of order 
of ahotC s as in an ignition front (Meyer 1984). We find no 
significant change in the speed of heating fronts at a given 
ftront when the mass of the central white dwarf is decreased 
from Mi = 1.2 M e to 0.6 M© in model W.l.cO.02. 

It was argued by Lin et al. (1985) that the velocity of a 
transition front must be close to the speed of gas in the front. 
We compared the speed of heating fronts with the maximum 
speed of gas inside the fronts, V r , ma x (see Figs, [l] and||). The 
velocity of inside-out and outside-in heating fronts appears, 
indeed, close to K,max in the fronts. For instance, Vfront ~ 
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3/2 V r ,max in models h0.1.c0.02 and 50hrl.5. Note that the 
speed of the cold gas into which a heating front penetrates 
is negligibly small compared to the speed of the front. 



3.4 Width of heating fronts 

The width of a transition front was defined in HMDLH as 
the region of the disc over which the variation of the viscosity 
parameter a = a(T c ) is 90 % of the total variation between 
Qcoid and «hot- This definition, although appropriate for an 
Qhot — Qcoid prescription, cannot be used when a is a smooth 
function of H/R (in particular model 50hrl.5). 

We instead define the transition region as a zone with 
large thermal imbalance where the temperature changes on 
short timescales. More specifically, if the local timescale of 
change in temperature is 



MR) 



d\nT c 



(R) 



(7) 



we define the transition region as the region of the disc where 
tT{R) < 3tT(i?front)- -Rfront is the location of the front, as 
defined previously (T c (i?f ront ) = 18,000 K). The factor 3 is 
arbitrary but it has been chosen to give results close to what 
an eye estimate would give. We confirmed that the above 
definition leads to widths which are comparable to those 
found in HMDLH but also comparable to those found with 
a definition of the transition region based on a deviation of 
the ratio \Q~ /Q + \ from unity. (This last definition should 
be equivalent to the definition used by, e.g., Cannizzo 1994; 
Cannizzo et al. 1995.) 

Figure |i| shows the fractional width (Sw / R{ Iont ) of heat- 
ing and cooling fronts (dots) as a function of the front loca- 
tion in the four models with various a-prescriptions. In each 
panel, the upper line of dots corresponds to cooling fronts 
and the lower line of dots corresponds to heating fronts. The 
width of cooling fronts is further discussed in §4.3. 

The width of a heating front is roughly proportional 
to the local vertical scale height H (note that H = c b /Q.k 
is defined at a fixed value T c = 1.8 x 10 4 K here). This 
is not unexpected since, except for the narrow precursor 
region, the thermal structure of the front is dominated by 
the heating term, Q + . Consequently, assuming that Vfront 
varies much less rapidly than iif^L (in agreement with our 
results, §3.3), Sw is the distance traveled by the heating front 
in a few thermal timescales Tth = (o^k) -1 (e.g. Frank, King 
& Raine 1992), i.e. 



Sw OC Vfront T t h OC (ac s /att K ) OC H OC -Rf/ont- 



(8) 



These results are consistent with the conclusions of Meyer 
(1984, 1986) and Papaloizou & Pringle (1985). 

A comparison of the various panels in Fig. ^ shows the 
influence of a on the width of a heating front. Increasing 
Qhot or decreasing a co id results in narrower fronts. It is not 
possible to interpret this result simply from Eq. (g), however, 
because the speed of heating fronts is also affected in a non- 
trivial way by changes in the viscosity parameter a (cf §3.3). 

In model 50hrl.5, transition fronts have profiles of 
8w/Rfiout which are slightly flatter than the solid line 
(H/Rf Iont ). According to the three other models, reduced 
widths are expected for larger values of a no t- The flatter 
profiles of Sw are therefore consistent, at least qualitatively, 
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Figure 4. The four panels show the fractional width (5w / Rf Tcm x) 
of cooling fronts and heating fronts in our four models with var- 
ious a-prescriptions (dots). In each panel, the upper line of dots 
corresponds to the width of successive cooling fronts and the 
lower line of dots corresponds to the width of successive heat- 
ing fronts. The solid line shows the fractional scale height of the 
disc {H/R{ Iont ). Note that heating fronts in model 50hrl.5 are of 
the outside-in type. 



with the additional variation of a with radius in this partic- 
ular model (q oc R 3 ' 4 ). 

We also investigated the effect of varying the mass of 
the central object on the width of transition fronts. We de- 
creased Mi from 1.2 M to 0.6 M in model h0.1.c0.02, 
which results in a reduced vertical gravity and therefore 
an increased vertical scale height H of the disc. At a given 
Vfront , the front widths are increased, but the ratio of heat- 
ing front (or equivalently cooling front) width to H is rea- 
sonably constant. This is additional evidence for a propor- 
tionality between Sw and the local thermal timescale in the 
disc, which depends on M\ via Q.k- 



3.5 Pairs of heating fronts 

When the thermal- viscous instability is triggered at a given 
ignition radius, two heating fronts appear which propagate 
in opposite directions in the disc. Because the ignition radius 
is often close to one of the disc edges, one of the two fronts is 
actually short-lived (i.e. quickly reaches the disc edge). The 
other heating front is long lived, propagates over a large 
range of radii and determines the global evolution of the 
disc. 

It is possible, however, to obtain models in which a pair 
of heating fronts appears somewhere in the middle of the 
disc. Figure [| shows such a situation (model 50hrl.5), where 
an inside-out and an outside heating fronts appear around 
R ~ 1.5 x 10 10 cm. In that case, the two fronts contribute 
significantly to the global evolution of the disc. We referred 
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Figure 5. Radial profiles of E and T c in a situation of coexistence 
of an inside-out and an outsidc-in heating front in the disc (model 
50hrl.5). Both fronts will significantly contribute to the global 
evolution of the disc. 

to this type of intermediate outburst as an outside-in out- 
burst. 

More generally, simulations show that the value of the 
ignition radius in a specific model depends on the mass 
transfer accretion rate Mt (HMDLH). This implies that 
a continuum of lightcurve shapes is possible, from quasi- 
symmetric (firm inside-out outburst) to strongly asymmet- 
ric (firm outside-in outburst). Consequently, the classifica- 
tion in two types of fronts initially proposed by Smak (1984) 
must be used with caution. 

A careful look at Figure 3 of HMDLH reveals a sudden 
increase and then decrease of the inner accretion rate, M acc , 
which precedes the bulk of accretion in inside-out outbursts; 
this spike is very short lived, contains almost no energy, 
and is probably not observable. This is a signature of the 
arrival of an in-propagating heating front at the disc inner 
edge. This front is the short-lived companion of the long 
lived inside-out heating front which propagates over a large 
range of radii in the disc. Outside-in outbursts (e.g. Fig. 4 
of HMDLH) do not show this burst of M acc because the 
outside-in heating front which reaches the disc inner edge 
is responsible for the bulk of the outburst accretion in that 
case. 

3.6 Reflection of an inside-out heating front 

The surface density after the passage of a heating front must 
be larger than the local value of E m i n (the minimum den- 
sity on the hot branch of the S-curve) to allow the front to 
propagate further in the disc. If not, a cooling front appears 
where E < E m i n in the disc (just behind the E spike) . Once 
the cooling front develops, the transport of excess angular 
momentum from the inner disc is strongly reduced and the 



Figure 6. Structure of a cooling front located at -Rf ro .nt ~ 3 x 
10 9 cm (model hO.l.cO.02). The structure is dominated by an 
outflow of gas and by radiative cooling. 

propagation of the heating front is stopped. The spike in 
E then diffuses under the action of viscous processes. This 
process is usually referred to as a reflection, although it cor- 
responds more exactly to the appearance of a new cooling. 

In the case of outside-in fronts, matter flows from the 
outer parts of the disc, which contain most of the disc mass, 
so that E steadily increases at any given radius, and the 
heating front is always able to reach the disc inner edge 
without being reflected in our models. 

The situation is quite different for inside-out heating 
fronts, for which mass is transferred from the inner parts to 
the outer parts of the disc. The comparatively small amount 
of mass in the inner disc leads to several reflections of inside- 
out heating fronts during one outburst cycle. After these 
multiple reflections, an inside-out heating front can reach 
.Rout without being reflected because enough mass has ac- 
cumulated in the outer disc and E behind the front remains 
above E m i n during the entire propagation. 

The reflections also occur in mini-outbursts (as com- 
pared to "normal" outbursts in which most of the disc - if 
not all - undergo a transition to the hot state) which are 
found in many numerical models (see e.g. Cannizzo 1993a, 
HMDLH). The significance of these mini-outbursts is un- 
clear because they have probably not been observed. 



4 COOLING FRONTS 

4.1 Structure of cooling fronts 

Figure ^| shows the structure of a cooling front located at 
fifront ~ 3 x 10 9 cm in model h0.1.c0.02. The innermost, 
hot, regions of the disc accrete at a much higher rate than 
the outer cold regions (see e.g. Fig. 9 of Cannizzo et al. 1995). 
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A substantial amount of angular momentum flows outward 
in the inner regions of the disk. The gas which carries this 
excess angular momentum is stopped when it reaches colder 
regions of the disc, where the transport of angular momen- 
tum is reduced. This results in a large outflow and an abrupt 
increase of E where rapid cooling sets in, as seen from Fig. ^| 

As noted by Vishniac & Wheeler (1996) (see also Vish- 
niac 1997), the inner disc is in quasi-steady state during the 
propagation of a cooling front; the accretion rate is almost 
independent of radius and it slowly decreases with time. The 
profiles of E and T c in these regions (not shown in Fig. ^[ 
but see, e.g., Fig. 5a of HMDLH) are close to those of the 
steady solution of Shakura & Sunyaev (1973). Between the 
inner quasi-steady disc and the front itself, there is a rela- 
tively broad region (called precursor region by Vishniac & 
Wheeler 1996) where the profiles of E and T c deviate sub- 
stantially from the Shakura-Sunyaev profiles. 

Figure |^ shows that the thermal structure of a cool- 
ing front is dominated, although not strongly, by radiative 
cooling (Q~): in first approximation, gas freely cools inside 
a cooling front (Vishniac & Wheeler 1996). The three non- 
local terms in the energy equation, namely Q a dv, Qj and 
Qpdv, are small but not negligible. The radial transport of 
energy (Qj) occurs in the region of strong gradient of T c . 
Advection of energy and adiabatic compression of the gas 
occur mainly in a region corresponding to the outflow of 
gas, i.e. where the gas which carries a substantial amount of 
heat (Q a dv) is compressed as it enters colder regions of the 
disc (Qpdv). 



4.2 Cooling front velocity 

Figure ^ shows the speed of cooling fronts in the four models 
with various a-prescriptions. Successive cooling fronts in an 
outburst cycle appear at various locations in the disc in mod- 
els using an a no t — a co id prescription. All the cooling fronts 
appear at the same location, close to the disc's outer edge, in 
model 50hrl.5. The initial speed of all cooling fronts is rel- 
atively high but soon, as they propagate inward, they relax 
to a reduced and common asymptotic speed which depends 
only on the current location of the front in the disc, Rfront, 
in agreement with early calculations by Cannizzo (1994). 

The effect of varying a co id on the (asymptotic) speed 
of a cooling front is small (at most). Increasing the value of 
Qhot results, however, in higher asymptotic speeds. 

In the asymptotic regime, Vfront has a weak depen- 
dence on -Rfront- This is especially true in models with an 
Qhot — «coid prescription (see also Cannizzo 1994). The dif- 
ferent variation of Vfront with i?t rorit in model 50hrl.5 can 
be accounted for by the variation of a with radius in this 
specific model, according to the scaling of Vf ront versus a no t 
in the three other models. Note that the sharp increase of 
Vfront close to the disc inner edge observed in all models is 
a boundary condition effect: the surface density E goes to 
zero at Ri n and this causes the front to propagate faster as 
E is closer to E m i n in the innermost regions of the disc. 

The speed of a cooling front, at a given radius R{ r 

ont i IS 

slightly sensitive to the mass of the central object. We ob- 
tained an increase of Vfront ~ 30% in the asymptotic regime 
when the mass M\ was changed from 1.2 Mq to 0.6 M©. 

Vfront in the asymptotic regime is not much larger than 
the speed of the hot gas into which the front penetrates; both 
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Figure 7. The four panels show the speed of successive cool- 
ing fronts in our four models with various a-prescriptions (dots). 
Cooling fronts appearing at various radii in the outer disc quickly 
converge to the same asymptotic speed (in a given model) as 
they propagate inward. (Note that all cooling fronts appear at 
the same location in model 50hrl.5.) The asymptotic speed of a 
cooling front slightly depends on the value of a n ot- 



are of order of a fraction of 1 km s _ . This is consistent with 
the idea that an evolution of the inner hot disc on a viscous 
timescale is typically required to allow the propagation of 
the cooling front (Lin et al. 1985, Vishniac & Wheeler 1996). 
But contrary to the case of heating fronts, we find no good 
correlation between the value of V r ,max inside the transition 
region and the speed of cooling fronts. 

Recently, Bobinger et al. (1997) directly determined the 
speed of a cooling front with the eclipse mapping technique. 
They observed the evolution of the brightness temperature 
profile in the disc of the dwarf nova IP Peg during four suc- 
cessive days of its decay from outburst. They infer a mean 
value of Vfront ~ 0.8 km s" 1 over the four days, which is 
consistent with the speeds shown in Fig. It will be diffi- 
cult, however, to determine the viscosity parameter a from 
such observations, even if variations of Vf r0 nt with Rf ron t 
were detected. For instance, models hO. l.OcO. 02 and 50hrl.5 
produce very different outburst cycles which could be dis- 
tinguished from observed lightcurves; yet the variation of 
Vfront with Rfront is not drastically different. A strong pre- 
diction of the DIM (which probably holds for any reasonable 
a-prescription) is, however, the rapid decrease of Vfront with 
Rfront shortly after the appearance of a cooling front. The 
detection (or non-detection) of this characteristic would be 
strong evidence (or disproof) that the DIM - as we under- 
stand it - operates in accretion discs. 

Vishniac & Wheeler (1996) proposed an analytical 
model of cooling fronts which uses the numerical results 
of Cannizzo et al. (1995) as calibration. In the derivation 
of Vishniac & Wheeler (1996), the speed of a cooling front 
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does not depend on the physical conditions in the cold gas 
behind the front. This is in agreement with the independence 
of Vfront with a co id observed in our simulations. 

Vishniac & Wheeler (1996) predict a speed for the cool- 
ing front 

, it, 0.54 

Vfront = OFCF i^J (9) 

if a = 50(H/ R) 15 (and assuming a Kramer's law opacity is 
valid for the hot gas), and a speed 



Vfront = "hot cf (~^r) 



(10) 



in the case of an a no t — "cold prescription. Here qf and 
cf are the viscosity parameter and the sound speed at the 
cooling front (i.e. where rapid cooling sets in). 

Equation (^) predicts a Vf ront which is slightly smaller 
(20% typically) than what we find in the asymptotic regime 
of model 50hrl.5 but the dependence with i?f ron t appears 
good. Equation (jTo|) predicts, however, a Vfront which is typ- 
ically twice as small as what is observed in the asymptotic 
regime of model h0.1.c0.02. 

Vishniac & Wheeler (1996) assume that the gas velocity 
at the cooling front, Vf, is ~ 1/6 olfCf and that Vf is much 
larger than the cooling front velocity. In our simulations, 
Vf ~ 1/7 cvfCf, but the cooling front velocity is not much 
smaller than Vf ■ The ratio of the gas velocity at the cooling 
front to the cooling front velocity is typically of order 2 (see 
Figs. [] and [?]). Neglecting Vfront is probably one of the reason 
for the discrepancy between the predictions of Vishniac & 
Wheeler (1996) and our results. 



4.3 Width of cooling fronts 

Figure i] shows the width of cooling fronts in the four mod- 
els with various a-prescriptions. Cooling fronts are much 
broader than heating fronts but their width Sw remains pro- 
portional to H . As for heating fronts, we interpret the width 
of a cooling front as the distance traveled by the front during 
a few thermal timescales (the thermal structure of the front 
being dominated by radiative cooling): Sw is proportional 



to Rli 2 nt cc H (Eq. (§)).[] 



This conclusion is in contradiction with the results of 
Cannizzo et al. (1995) who find Sw = V HR in a numerical 
model using a = 50(H/ R) 15 . Theprofiles of Sw/Rf Iont as 
a function of Rhoat shown in Fig. ^ for model 50hrl.5 are 

* Naively, one would expect from Eq. (^) smaller Sw for cooling 
fronts than for heating fronts (because of smaller Vf rimt ). Equa- 
tion (^) does not take into account, however, several complications 
in the problem. The thermal timescale r t h = (oOjf)" 1 ! which is 
very different for annuli in the hot and the cold regions of the 
disc, is probably increased because of the ionization energy in 
transition regions (Papaloizou & Pringle 1985). The precursor re- 
gion drives out of thermal equilibrium the disc annuli entering the 
transition region in the case of heating fronts (see Fig. ^), while 
there is no equivalent for cooling fronts. Finally, radiative cooling 
clearly does not dominate the thermal structure of a cooling front 
as much as viscous heating does in a heating front, which proba- 
bly leads to comparatively longer thermal timescales for cooling 
fronts. A key feature that Eq. feh includes, however, is the strong 
variation of the thermal timescale with radius. 



Si 





Figure 8. Successive normalized profiles of S and local mass 
flow rate M as a function of the dimensionless radius (R — 
R'm)/(Rfront ~ Rin) during the propagation of a cooling front 
(model hO.l.cO.02). The superposition of profiles reveals the 
nearly self-similar evolution of the inner hot disc. 



consistent with those of Cannizzo et al. (1995). Cannizzo 
et al. (1995) claim that Sw = \J HR is also true in models 
using an dhot — « C oid prescription. Our three models with 
Qhot — Qf co id prescriptions rule out this possibility and show 
that Sw oc H. As mentioned previously, we interpret the 
flatter profiles of Sw in model 50hrl.5 as due to the varia- 
tion of a with radius in this particular model. Models using 
a = oto(H/R) n with various n produce various slopes for 
the profiles of Sw as a function of R{ Ton t ■ 

A comparison of the various panels in Fig. ^| shows that 
the effect of varying a on the width of a cooling front is 
small. 



4.4 Self-similar solutions 

Vishniac (1997) proposed a self-similar solution for the 
structure of the hot inner disc during the propagation of 
a cooling front. We confirm the existence of a self-similar 
regime, although the self-similarity described here is some- 
what different from the solution of Vishniac (1997) because 
E is found to scale naturally with S m i n . We leave a detailed 
analysis of the self-similar solution for a future paper, but 
we give here our main qualitative results. 

Figure ^ shows successive profiles of E and local mass 
transfer rate M during the propagation of a cooling front in 
model W.l.c0.02, when the asymptotic regime is reached. 
The profiles of E have been normalized to the minimum 
surface density in the cooling front, which is extremely close 
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to E m in at -Rfroaf This is not unexpected, since the transi- 
tion from the hot to the cold state occurs precisely when E 
reaches E m i n . The mass flow rate M has been normalized, 



following Vishniac (1997), to the mass transfer rate at the 
disc inner edge, M acc . The horizontal coordinate chosen here 
is (R — Rin) / (Rhout — Rin) in order to minimize the effect 
of the inner boundary condition on the profiles. The value 
of flfront varies by more than an order of magnitude for the 



8l and yet the normalized E 
ar: except for the effect of 



successive profiles shown in Fig. 
and M remain remarkably simi 
boundary conditions, the evolution of the inner hot disc is 
close to self-similar. 

The mass accretion rate onto the compact object, M acc , 
is proportional to Mf ront at the transition radius. But 



Mfront = 27Tiif ron tE(J?f r ont)K-(-Rf] 



(11) 



ont ) — S m in (Bfat) oc R^ nt (e.g. HMDLH), and 
Vr(-Rfront) is comparable to the front velocity. If a in the 
hot inner disc is constant, one expects a weak dependence 
of K(-Rfront) with .Rfront and M acc oc R" Iont , with n slightly 
larger than 2. Indeed, numerical simulations show that 
n ~ 2.2 for an ahot — Q co id prescription, and n ~ 2.4 when 
a oc (H/R) 15 . 

As discussed by Vishniac (1997), the hot inner disc 



evolves on a viscous time during the propagation of a cool- 
ing front, and is therefore very close to thermal equilibrium. 
The inner disc empties as a result of accretion onto the com- 
pact object, but also because an even larger amount of mass 
is transferred into the cold outer regions of the disc. These 
parts of the disc are essentially frozen on the time scale of 
propagation of the cooling front. Consequently, E(i?) behind 
the front remains constant in time and is equal to if E m i n (i?) , 
where if is a self-similarity constant that depends on the 
mass of the primary and the viscosity parameter a both in 
the hot and the cold regions of the disc. For the parameters 
of Fig. |^, K ~ 4. The constant K increases with the mass 
of the primary (for example, K ~ 6 — 7 for a 7 Mq primary 
typical of a black hole SXT). 

The sudden increase of E just after the passage of a cool- 
ing front is actually smaller in the pre-self-similar regime, so 
that E/E m i n increases as the front propagates toward the 
asymptotic regime. This accounts for the density profile in 
the outer disc which is relatively flat close to the outer edge, 
steepens at smaller radii and finally becomes proportional 



to E n 



in the self-similar regime. 



4.5 Reflection of a cooling front 

The sudden increase of E just after the passage of a cool- 
ing front cannot be arbitrarily large (Vishniac, 1997). If 
E > E max behind the cooling front, an inside-out heating 
front develops and propagates over some distance in the disc 
because E is substantially larger than E m i n in the outer disc. 
The initial cooling front soon disappears as hot gas from the 
outer regions of the disc accretes inward. 

This situation corresponds to the reflection of a cooling 
front into a heating front and it occurs whenever E max /E m i n 
is too small. The ratio E max /E m j n has a weak radial depen- 
dence (#0.02-0.03. e g HMDLH). If the self-similarity con- 
stant K is too large, the cooling front will experience a re- 
flection typically when it enters the self-similar regime: the 
gradual increase of E/E m i n preceding the self-similar regime 



will drive E behind the front above E max and lead to the re- 
flection. 

The multip le reflections of cooling fronts are well known 
if a co id = Qhot (Smak, 1984). In this case, the inner disc can 
never be brought entirely into quiescence, and the reflected 
transition fronts propagate back and forth in a restricted 
region of the disc, causing small amplitude quasi-oscillations 
of the disc luminosity. 

Note that because K increases with the mass of the 
primary, the reflections of cooling fronts are more likely to 
occur in BH SXTs than in neutron star SXTs or dwarf novae. 



4.6 Exponential decays 

We reexamine, in light of our results on cooling fronts, the 
claim by Cannizzo et al. (1995; see also Vishniac & Wheeler 
1996) that a = ao(H/ R) is necessary to account for ex- 
ponential decays of soft X-ray transients (SXTs; but see also 
the similarity between SXTs and dwarf novae pointed out 
by Kuulkers, Howell & van Paradijs 1996) 

Our simulations show that the inner accretion rate 
M acc oc Rf Iont in the self-similar regime of propagation of a 
cooling front (n = 2 — 2.5 depending on the a- prescription). 

Since X-rays are produced close to the central object 
in an SXT, the exponentially decaying lightcurves imply an 
exponential decay of M acc as well, i.e. 



dM a , 



dt 



OC M a , 



(12) 



If we assume that the decay phase corresponds to the prop- 
agation of a cooling front in the disc and that our results 
apply (no irradiation effect is included in our calculations), 
then 



dM a , 



dt 



nll-l T/ n> Mfront 

OC tt bont Vfront OC 



(13) 



and an exponential decay is possible only if Vf ront /Rkont 
remains constant, meaning that ifeont must also vary expo- 
nentially with time: 



Vfront — 



dRhon 

d~r 



OC iifron 



(14) 



Note that this requirement was already deduced, although 
differently, by Cannizzo et al. (1995) and Vishniac & 
Wheeler (1996). 

Figure ^ shows that Vf ron t is primarily independent of 
iifront (in the asymptotic regime) in the three models with 
Qhot — «coid prescriptions. On the contrary, model 50hrl.5 
shows a variation of Vf ront with i?f ront which is linear in 
first approximation. Our calculations therefore confirm that 
a dependence of the viscosity parameter a with radius (close 
to a oc R 3 / 4 like in model 50hrl.5) is required to obtain 
exponential decays of M acc in the DIM. 

A major concern however is the relevance of interpret- 
ing exponentially decaying lightcurves as due to the propa- 
gation of a cooling front in the discs of SXTs. In particular, 
Shahbaz et al. (1998) argue that both linear and exponen- 
tial decays are observed in SXTs (see also the compilation 
of lightcurves by Chen et al. 1997). Since there is no obvious 
physical reason why a would vary markedly from one system 
to another, the possibility that both linear and exponential 
decays are observed in SXTs argues against an interpreta- 
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tion of the exponential decays as due to a specific functional 
form of a. 

An alternative and perhaps more promising explanation 
for the exponential decays has been proposed by King & 
Ritter (1998). The authors argue that strong irradiation of 
the disc may prolongate the phase during which the entire 
disc is in the hot state and lead to an evolution of the disc 
on a (long) viscous time (see however Dubus et al. 1998). 

4.7 Cooling fronts and the periodicity of outburst 
cycles 

The asymptotic regime of cooling front propagation has a 
profound influence on the outburst cycles experienced by the 
disc. Since the inner hot disc evolves in a self-similar way, the 
profiles of E and T c inside -Rf ront are uniquely determined 
by the current value of -Rf ron t • 

A cooling front following a heating front that reached 
the outermost regions of the disc therefore propagates in the 
exact same way a previous cooling front starting from Rout 
did. This type of cooling front erases the history of the disc 
and is responsible for the periodicity of the outburst cycles 
in the simulations. An outburst cycle ends and a new cycle 
begins when a heating front reaches R ou t (also corresponding 
to the largest amplitude outbursts). 

Observed outburst cycles in dwarf novae are usually reg- 
ular but not periodic (e.g. Warner 1995, Cannizzo & Mattei 
1992). This may hint that a piece of physics is missing in 
the standard DIM. A variation of the mass transfer accre- 
tion rate Mt appears as a natural candidate for explaining 
the irregularity of observed cycles. 



5 EVOLUTION AROUND S-CURVES 

Figure |^ shows limit cycles experienced by two disc annuli in 
models h0.1.c0.02 and 50hrl.5 in the E — T e s plane, where 
T e ff is the effective temperature of the disc. 

Global effects (i.e influence by neighboring annuli) ap- 
pear when an annulus no longer lies on the S-curve. The 
passage of a cooling front corresponds to a jump from the 
upper hot branch of the S-curve to the lower cold branch. 
Inversely, the passage of a heating front corresponds to the 
jump from the lower branch to the upper branch. Note that 
the maximum value of E reached by a disc annulus on the 
lower cold branch of the S-curve during a limit cycle is in- 
dicative of the proximity of this annulus to the ignition ra- 
dius in the disc. 

Figure ^ shows the evolution of disc annuli during the 
passage of a transition front: a sudden increase of E is as 
good a signature of the transition as a change of T e g is. For 
a cooling front, the sudden change of E shortly precedes 
the phase of rapid cooling of the annulus. For a heating 
front, the sudden change of E corresponds to the passage 
of the spike of E. This shows that inside a transition front, 
the viscous time is not much larger but comparable to the 
thermal time since surface density and effective temperature 
vary on similar timescales. 

The complexity of the limit cycles experienced by the 
annulus at R — 10 10 cm in model h0.1.c0.02 is due to 
the many reflections of inside-out heating fronts into cool- 
ing fronts around 10 10 cm in this model. Note that the 
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Figure 9. Limit cycles experienced by disc annuli located at 
R = 10 9 and 10 10 cm in the £ - T off plane. The left pan- 
els correspond to model h0.1.c0.02 and the right ones to model 
50hrl.5. The solid lines are the thermal equilibrium curves, where 
Q+ = Q~ (standard S-curves). The dashed lines trace the suc- 
cessive states of a disc annulus, while the density of circles is 
representative of the amount of time spent by an annulus in a 
given state. The evolution of the disc annulus shown in the lower 
left panel is complex because this annulus is located in a region 
of front reflections. 

track followed by the annulus states during the successive 
jumps from the upper branch to the lower branch do not 
cover each other in that case because cooling fronts have 
not yet reached the asymptotic regime at 10 10 cm in model 
hO.l.cO.02. 

Note that the evolutions around S-curves shown here 
differ significantly from the results of Ludwig, Meyer- 
Hofmeister & Ritter (1994) using the theory of infinitely 
thin transition fronts (Meyer 1984). Some additional limi- 
tations of the theory were pointed out by Ludwig & Meyer 
(1998). 



6 VALIDITY OF THE DISC EQUATIONS FOR 
THE TRANSITION FRONTS 

It is standard in the theory of thin accretion discs to assume 
that the vertical structure of the disc is decoupled from its 
radial structure because the former evolves on much shorter 
timescales than the latter (e.g. Frank et al. 1992). Because 
E and T e g vary on similar timescales, this assumption is no 
longer valid inside transition fronts: the vertical structure is 
coupled to the radial structure of the disc. In particular, ra- 
dial transport of energy becomes important in the transition 
regions, while the transport occurs only vertically in steady 
thin discs (Frank et al. 1992). 

A correct treatment of the coupling between the radial 
and vertical structures of the disc would probably involve 
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Figure 10. An illustration of the sensitivity of the propagation 
velocity of heating fronts to the poorly known efficiency of radial 
transport of energy (term J in Eq. (2)), in model hO. I.c0.02. 

2D simulations. For simplicity, however, the radial trans- 
port of energy is usually parameterized in a simple manner 
in ID disc equations (see also the approximations of decou- 
pling used for the determination of the cooling rate Q~ , 
HMDLH). The fact that various prescriptions exist for J 
(Eq. (^)) shows that the radial transport of energy in the 
disc is not well known, and surely not within a factor of a 
few. 

Figure [To| shows the effect on the speed of heating fronts 
of increasing or decreasing by a factor 3 the intensity of 
the term J. The large uncertainty in J translates into a 
significant uncertainty on the speed of heating fronts. (Note 
that the speed of cooling fronts is not affected because the 
radial transport of energy has a minor influence on their 
structure; cf Fig. fy. 

An additional concern about the validity of the disc 
equations inside transition fronts comes from the narrowness 
of the fronts. Lin et al. (1985) pointed out that having a 
front width Sw ~ H, although it does not affect much the 
Keplerian character of the flow, could lead to a Rayleigh- 
unstable situation in the disc. A comparison between the 
criterion derived by Lin et al. (1985) and the structure of the 
narrow heating fronts found here suggests that this type of 
fronts may indeed experience the development of a Rayleigh 
instability. This could result in an enhanced radial transport 
(Lin et al. 1985) and an additional uncertainty in J. 

The equations used in this study are only valid in the 
thin disc approximation, i.e. when radial pressure gradients 
can be neglected and both and dQ/dR are equal to the 
Keplerian values. The effects of deviations from Keplerian- 
ity and radial pressure gradients were shown, for instance, 
by Ludwig & Meyer (1998) using full hydrodynamical disc 
equations and by Kley & Lin (1996) in boundary layer com- 
putations. Although these effects become non-negligible in 



narrow heating fronts, we do not expect our main results 
on the structure and properties of transition fronts to be 
affected if they were fully taken into account. This could be 
checked by taking the full equations, which is beyond the 
scope of the present paper. 

Finally, we note that the use of a local a-prescription in 
narrow transition regions (8w ~ H) is also subject to cau- 
tion because the approach of scale separation becomes much 
less rigorous on scales so close to the size of the "turbulent 
eddies" . 



7 DISCUSSION 

Our simulations show that transition fronts systematically 
experience a deceleration during their propagation in the 
disc, whatever their direction of propagation. Yet, the front 
widths increase with radius in all the cases. This shows that 
the speed of a front is not directly related to its width (see 
also a related claim by Vishniac & Wheeler 1996 for cool- 
ing fronts) . A description according to which the width of a 
transition front, which defines the gradients and the fluxes 
in the front, also determines its speed thus appears oversim- 
plified. 

Our simulations reveal a qualitative agreement with one 
of the predictions of the theory of infinitely thin transition 
fronts (Meyer 1984, 1986; see also Papaloizou & Pringle 
1985): we clearly observe that the speed of transition fronts 
is related to the proximity of E in the disc to the critical val- 
ues E m i n (for cooling fronts) and E max (for heating fronts). 

Although Ludwig et al. (1994) pointed out that this 
theory may not be applicable to relatively broad cooling 
fronts, it could provide an explanation for the rapid decel- 
eration of cooling fronts: the values of E around the region 
where the cooling front appears are close to E m i n , so that 
the front speed is initially high; soon, however, the front en- 
ters regions with larger and larger values of E as compared 
to E m i n (quasi-Shakura-Sunyaev profile), which results in a 
gradual decrease of the speed; later, the front must wait for 
the inner hot disc to evolve on a viscous timescale, typically, 
before E approaches E m i n in the disc, which explains the 
small speeds in the asymptotic regime. On the contrary, the 
speed of a heating front is able to remain large during the 
propagation, possibly because the profile of E encountered 
in the disc is typically in between E m i n and E max (Cannizzo 
1993b), i.e. not too far from E max everywhere in the disc. 
If this is true, the deceleration of heating fronts is due to 
the (necessary) proximity of E to E max close to the ignition 
radius but not elsewhere in the disc. 

We have not considered the structure and properties of 
transition fronts close to Rout in this paper. Indeed, addi- 
tional simulations show that in a small region close to the 
disc outer edge, transition fronts are affected in a complex 
way if Rout is allowed to vary with time. The behavior of the 
fronts in this region of the disc may also depend sensitively 
on the precise choice of the outer boundary condition and 
how mass is deposited in the disc. 
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8 CONCLUSION 

In this paper, we have investigated the structure and proper- 
ties of transition fronts in thin accretion discs with detailed 
numerical calculations. 

We showed that heating fronts are very narrow and have 
complex structures. They propagate at a speed which de- 
pends on the profile of surface density in the disc, the radial 
transport term in the energy equation and the value of the 
viscosity parameter a, but is typically of order a few km s _1 . 

Cooling fronts are broader, have a simpler structure and 
have smaller speeds (of order a fraction of a km s _1 ) than 
heating fronts. We found that their width is not equal to 
V HR, but is rather oc H, the local disc scale height, as for 
heating fronts. 

We confirmed that the structure of the inner hot disc is 
well described by a self-similar solution during the propaga- 
tion of a cooling front. We proposed such a solution in which 
the surface density does not scale arbitrarily but with the 
value of the critical surface density E m i n . The self-similarity 
of the disc appears responsible for the periodicity of the out- 
burst cycles in the simulations. 

Since all our models show a deceleration of transition 
fronts during their propagation in the disc, the observation 
of such a deceleration would constitute a nice confirmation 
that the thermal-viscous disc instability, as we understand 
it, is responsible for the large amplitude variability of discs 
around white dwarfs, neutron stars and black holes. 
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